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ABSTRACT 

We analyze the electromotive force (EMF) terms and basic assumptions of the linear and nonlinear 
dynamo theories in our three-dimensional (3D) numerical model of the Parker instability with cosmic 
rays and shear in a galactic disk. We also apply the well known prescriptions of the EMF obtained by 
the nonlinear dynamo theory (Blackman & Field 2002 and Kleeorin et al. 2003) to check if the EMF 
reconstructed from their prescriptions corresponds to the EMF obtained directly from our numerical 
models. We show that our modeled EMF is fully nonlinear and it is not possible to apply any of 
the considered nonlinear dynamo approximations due to the fact that the conditions for the scale 
separation are not fulfilled. 

Subject headings: ISM: galactic dynamo — magnetic field 



1. INTRODUCTION 

It seems that the issue of the magnetic field am- 
plification in galaxies may be well explained by the 
two main physical mechanisms: the Parker instability 
(PI), which takes into account the cosmic rays (CR) 
and t he shear (e.g. lHanasz fc Lesch 1 12003: Ha nasz et al. I 
12004 and references therein), and the magneto- 
rotat i onal instability (MRI, e.g. i Dziourkevitch et al. I 
12004 iKitchatinov fc Riidiger ll2004f ). The possible sce- 
nario of the magnetic field evolution could be presented 
as follows: when the protogalaxy starts to rotate differ- 
entially, the MRI mechanism occurs, and this results in 
a very efficient magnetic field amplification even to the 
level of ii G with the global e-folding t ime of 100 Myr or 
even less (Dziourk evitch et al. 1 12004 ). Simultaneously, 
the quadrupole symmetry of the large-scale magnetic 
field is being created. MRI also causes the turbulent mo- 
tions in the galactic disks. The process of supernovae 
(SN) explosions, th at arises in young galactic objects 
(e.g. IWidrow 1 120021 and references therein) suppresses 
the MRI mechanism. In the same time, the cosmic rays 
produced in SN re mnants may induce the Parker insta - 
bility process (e.g. iParker I Il99l lHanasz et aT7~l l2004f ). 
Hence, we may conclude that during the early stage of 
galaxy evolution the MRI process is being replaced by 
the PI mechanism. 

The local simulations of the large-scale magnetic field 
took into account the turbulent dynamo theory and 
the magnetic back- r eaction onto the turbulent motions 
dPiddington I [liM )l97l fl97fl iKulsrud fc Anderson I 
Il995f h The authors drew the conclusion that it was 
difficult to obtain the am plification of the total magnetic 
field (e.g. IWidrow 1 l2002t and references therein). It 
might be explained by the equipartition of the random 
magnetic field component with the random turbulent 
motions. Such process suppresses the dynamo action at 
later times. Moreover, the magnetic field amplification 
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might be easily stopped even by the presence of the 
weak l arge-scale magnetic field I Cattaneo fc Vainshtein I 
1991t IVainshtein fc Cattaneo I Il992t ICattaneo I 11994 . 
Cattaneo fc Hughes I [l996t IZiegler et al~l fl99ll . 
iBlackman fc Fieldl (|1999t ) explained that papers 
analyzing the analytically strong suppression of the 
dynamo coefficient a should distinguish between differ- 
ent state orders of the turbulent quantities. However, 
their analysis did not completely solve the problem 
of the quen ching of the dynamo coe fficients. In their 
next paper ([Blackman fc Field 1120001). they prov ed that 
the results from the ICattaneo fc Hughes I (l996) model 
were based on the assumption about the periodicity of 
the boundary conditions. Nevertheless, the quenching 
effects could also appear even when the open boundary 
cond itions were applied (e.g. iBrandenburg fc Dobler I 
12001 . 

The classical dynamo theory doe s not conserve the 
total magnetic helicity (see e.g. [Blackman fc Fieldl 
120021 BF02). In media characterized by the high 
magnetic Reynolds number (R m 3> 1) the total 
helici ty should remain c o nstant in closed regions 
e.g. iBerger fc Field I 11984 IBrandenburg et ai~l 120021: 



Brandenburg fc Subramanian I 120051 : 



III. 



ISubramanian I 

The permanent helicity is also an ad- 
ditional factor , which suppress es th e dynamo 
activity (e.g. IBrandenburg et aT~l I2002T). The 

following p apers (IBlackman fc Field I 120001: 

Kleeorin et all 120001: 



Blackman & Brandenburg 



2003 1 : IKleeorin fc Rogachevskii 1999t IKleeorin et al. 



2002L I2000L 120031 : iRogachevskii fc Kleeorin 1 120001 l200l 



presented the two methods that allowed the modeling of 
the dynamo action evading the problem of the constant 
helicity. The first method is based on the ejection of 
the magnetic helicity through boundaries. The second 
one uses the creation of the negative and positive 
helicity at the large and small scales respectively (se e 
Brandenburg et al. 1 12002) : IKleeorin et £d~f I2002L [2003). 



Blackman fc Field I ()2002f ) in their next paper analyzed 



the nonlinear prescription of both dynamo coefficient 
a and (3. Both factors were obtained without any 
linearization and took into account all terms in the 
equation of the evolution of the fluctuating part of the 
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magnetic field (see Eq. [7]) . The results were similar to 
the d ynamo coefficients obtained by (|Pouquet et al. I 
11976ft . but the units were different (without the time 
integration). They also solved the small- and large-scale 
helicity dynamo equations numerically simultaneously 
with the equation for the EMF time evolution. That 
allowed them to obtain the growth of the large-scale 
magnetic field in the kinematic phase. 

The new form of the dynamo coefficients for 
anisotropic turbulent motions with the presence of 
the large-scale magnet i c fiel d was presented by 
iRogachevskii fc Kleeorin I ((2001). They calculat ed dy- 
namo coefficients accor ding to the IRadler I (|1980ft EMF 



2005), which ne- 



prescription (see also iKowal et al. 
glecte d all quadratic t erms in the mean field in the 
EMF. IKleeorin et al. I (|2003f ) used those forms of the 
dynamo coefficients to solve numerically the dynamo 
equation in the local thin-disc approximation. The au- 
thors took into account the quenching of both coeffi- 
cients, a and /3, and helicity flux through the boundary. 
They found that it was possible to obtain the growth 
of the large-sc ale magnetic field whe n a-quenching was 
only analyzed ((Kleeorin et al. I f2002) . If the model in- 
cluded also /3-quenching, no growing solution of the dy- 
namo (|Kleeorin et al. 112003ft could be obtained. Both 
in Blackman-Field and in Kleeorin-Rogachevskii ap- 
proaches there is an a m term that quantifies the small 
scale helicity current. This term depends on the current 
helicity flux, and there are different theories for this flux. 
In Kleeorin et al. a heu ristically motivated expression 
for th e flu x was used, in Brand enburg fc Subramanian I 
(12005ft thelVishniac fc Cho I (|2001ft flux was used, and in 
iShukurov et al. I ( 2006ft a simple advective flux was used. 
In all these cases the current helicity flux allows the field 
to saturate at high levels. 

The latest research on the turbulent enforcement 
in the solar convective zone calculated in the lo- 
cal cube with the shear has shown that the open 
boundaries help to obtain the amplification of the 
large-scale magnetic field even without the helic- 
ity of turbulent motions dBrandenburg et al. I 120051 : 
IBrandenburg fc Subramanianl 120051 end references 
therein). We have to stress that their result, an increase 
of the total magnetic energy, was obtained without 
any assumption considering the additional EMF of 
dynamo. The authors applied isotropic and homoge- 
neous turbulence with and without the helical forcing 
in their model. On the other hand, Brand enburg and 
his collaborators (|Brandenburg et al71l2005ft interpreted 
their results in terms of the mean field dynamo the- 
ory. The time evolution of a was obtained from the 
calculated electromotive force (|Brandenburg fc Sandin I 
20041: IKowal et al. Il2005|: IBrandenburg fc Sokoloff 1120021 ; 



Brandenburg Tl2001ft~ IBrandenbur g fc Sandin I (2004 ) 
explained that thanks to the flux of the current helicity 
flowing out of the cube the process of a-quenching tends 
to be not as disastrous as forseen. The only conditions 
are the intermediate level of R rn and open boundaries. 
When the high value of R m is applied to the model, the 
comparatively lower value of the large-scale magnetic 
field strength are obtained (IBrandenburg fc Sandin I 
2004). However, in astrophysical objects R m is always 
high. That is why this result seems to be peculiar. On 
the other hand, the authors made it clear that the total 



magnetic energy grows mainly in the kinematic phase 
of the dynamo and their results do not depend strongly 
on R m . Futhermore, they calculated the value of a 
coefficient based on the modeled EMF. The obtained 
factor was similar to the same c oefficient cal c ulated 
according to the IKleeorin et al~l (|2000l . 120021 12003ft 
prescriptions. We believe that the fact that such simi- 
larity occurs results from the isotro pic and homogeneous 
turbu lence in both models (IBrandenburg et al. I 

2005: IBrandenburg fc Subramanian I 120051 : 

Brandenburg fc Sandin I l2004ft . It may also happen 



due to the fact that the authors applied standard 
dynamo approximations, which neglect the quadratic 
terms in the mean field in EMF. 

The previously mentioned results indicate that the 
realistic physical simu lations are of the great im por- 
tance wh en the MRI (fDziourkevitch et al 1 120041) and 
the PI (lHanasz et al. 1 120041: lHanasz fc Lesch I 120031: 
IKowal et al. 112005ft processes are considered. Both mod- 
els, which meet enumerated requirements, showed that it 
was possible to amplify galac tic magnetic fie l d efficiently . 
lHanasz fc Lesch I (|2003ft and lHanasz et al. I ((200412005ft 
presented that the following two processes: the Parker 
instability driven by cosmic rays from supernovae and 
the shear from the differential rotation enable the mag- 
netic field amplification (with the e-folding time scale of 
250 Myr or even 140 Myr). T he model also a pplied re- 
alistic gravity according to the iFerriere I ((1998ft prescrip- 
tions. 

The idea of obtaining the dynamo coefficients (a 
and f3) from the electromotive force calculated in 
the local numerical simulations pr oved to be es s ential 
for many other authors too ( e.g. IZiegler et al. 1 119961 : 
IBrandenburg fc Sokoloff] l2002t etc). We included the 
calculations of the dynamo coefficients, which we ob- 
tained from the calc ulated EMF in our previous paper 
((Kowal et al. 1 12005ft . The application of the statisti- 
cal methods provided us with the acceptable values of 
the dynamo a-tensor. On the other hand, the values of 
P coefficients were negative. Such values are in consis- 
tent with the Radler prescription ((Radler 1 11980ft . This 
may be caused either by the applied statistical method, 
which does not take into account physic al differentia- 
tion, or by the linear EMF approximation ((Kowal et al. I 
2005). For this reason we decided to analyze that mat- 
ter in our present study. We search for the condi- 
tions, which should be fullfilcd in order to make lin- 
earization of the electromotive force in the mean field 
dynamo theory possible. We would like to examine 
the following problems: the scale separation, the ra- 
tios of the terms in the equ ation for the small-scale 
magnetic field evolution (e.g. IRadler 1 11980ft . the mag- 
nitude of the turbulent kinetic energy in comparison 
to the large-scale magnetic one. Next, we plan to ap- 
ply ^_the_e^timatiojisof the dynamo coefficients derived 
bvlBlackman fc Field I (l2000ft: IRogachevskii fc Kleeorin I 
((2001ft and IBrandenburg fc Sandin I ((2004ft to our mod- 
els. We would like to check if their approximations fit 
into the calcu l ated electromotive force in our models 
lHanasz et al. I ((2004ft . In this part of this work we do 
not include the explicit analysis involving the conserva- 
tion of the magnetic helicity. The investigation of the 
magnetic helicity conservation in our models is already 
advanced, but it is complex enough to be described in 
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separate paper, which is unde r preparation. We can- 
not ap ply t he approximations o f lRuec iger fe Kichatinovl 
(1993) and iKitchatinov et al. I (j!994l ) derived from the 
EMF quenching by the usage of the Second Order Corre- 
lation Approximations. They assumed that the Strouhal 
number S is essentially smaller than 1 (S <C 1, where 
S = t c x v/l c ). This assumption is not fullfiled in our 
numerical experiments, where S is about 1. Finally, we 
discuss our results. 

2. NUMERICAL MODELS OF THE COSMIC-RAY DRIVEN 
DYNAMO 

The first complete 3D numerical mod el of the CR- 
driven dynam o has been demonstrated bv lHanasz et al.~l 
(|2004 120051) . The presented numerical simula- 
tions were performed w ith the Zeus-3D MHD code 
(iStone fc Norman I ll992allbT) extended with: the cosmic 
ravs (|Hanasz fc Lesch 1 12003). supplied by the supernova 
remnants randomly exploding in the disk volume, the 
resistivity of the ISM leading to the magnetic recon- 
nection, the shearing-periodic boundary conditions, the 
rotational pseudo-forces and the realistic vertical disk 
gravity. The full set of equations describing the model 
includes the set of resistive MHD equati ons completed 
by th e cosmic ray transport equation (see lHanasz et aT~l 
l2004f) . where anisot ropic diffusi o n of c osmic rays is im- 
plemented following iRvu et al~l (|2003[ L 

The principle of the action of the CR-driven dynamo 
is based on the cosmic ray energy supplied continuously 
by SN remnants. Due to the anisotropic diffusion of cos- 
mic rays along the horizontal magnetic field lines, cosmic 
rays tend to accumulate within the disc volume. How- 
ever, the configuration stratified by the vertical gravity 
is unstable against the Parker instability. Buoyancy ef- 
fects induce the vertical and horizontal motions of the 
fluid and the formation of the undulatory patterns - 
magnetic loops in the frozen-in, predominantly horizon- 
tal magnetic fields. The presence of rotation in galactic 
disks implies a coherent twisting of the loops by means 
of the Coriolis force, which leads to the generation of the 
small-scale radial magnetic field components. The next 
phase is merging the small-scale loops by the magnetic 
reconnection process to form the large scale radial mag- 
netic fields. Finally, the differential rotation stretches the 
radial magnetic field to amplify the large-scale azimuthal 
magnetic field. The coupling of amplification processes 
of the radial and azimuthal magnetic field components 
gives rise to an exponential growth of the large scale mag- 
netic field. The timescale of the magnetic field amplifi- 
cation, resulting from the action of CR-driven dynamo, 
has been foun d to be about 250 M yr in typical galac- 
tic conditions (jHanasz et al. Il2004h . close to the galactic 
rotation timescale. 

In our model of th e cosmic rays driven Parke r in- 
stability with shear (jHanasz et al. I 12004 l2005h we 
have included the following physical elements: the 
cosmic ray component descri bed by the diffusion- 
advection transport equation (see lHanasz fc Lesch 1 12003. 
for the details of the numerical algorithm), cosmic 
rays diffusing anisotropical l y along magnetic fi eld lines 
(jGiacalone fc Jokipii I IX9Q9I : I Jokipii et al. Ifl999h . super- 
nova r^nrnants_e2cpjpdingrandomly in the disk volume 
(see lHanasz fc Lesch 1 [2000) , the finite (currently uni- 
form) resistivity of the ISM (see lHanasz et al. 1 120021 : 



lHanasz fc~ Lesch 2003; Tanuma et al. I 



2003 ) and the re- 
alistic vertical disk gravity (jFerriere 1 11998 ). The sys- 
tem of coordinates x, y, z corresponds locally to the 
global galactic cylindrical system r,4>,z. The disk ro- 
tation was defined by the values of the angular velocity 
£1 =0.05 Myr -1 and the value of the shearing param- 
eter is q = 1. The boundary conditions are periodic 
in the Y-direction, s heared in the X-direction (following 
lHawlev et ai~lll995| ) and open in the Z-direction. 

The open boundary conditions in the Z-direction are 
constructed in such a way that fluid, along with the mag- 
netic field frozen in the fluid, can move out through the 
upper and lower box boundaries. These boundary condi- 
tions rely on the copying of fluid variables and magnetic 
field components from the last cell-layer inside the com- 
putational domain to the ghost zones. 

In this paper we analyze the numerical mode l of 
cosmic-ray driven dynamo by lHanasz et al. I (|2004l ) to- 
gether with a series of the new three models that form a 
subset of the new extensive series, w hich are described in 
detai l in two complementary papers (jHanasz et al. Il2006l . 
2007). In the presented simulations we adopt the basic 
parameters (the vertical gravity model and the thermal 
gas column density) from the global model of ISM in 
Milky Way (jFerriere II 19981) for the galact ocentric radius 
R=8.5 kpc (model A lHanasz et al. Il2004h and R=5 kpc 
(new models B, C and D). Moreover, models B, C and 
D rely on improved boundary conditions for cosmic rays 
(with respect to Model A) , imposing e cr = on outer Z- 
boundarics. The initial magnetic field in the azimuthal 
direction corresponds to a small fraction 10~ 8 of thermal 
energy for model A and 10 -4 for models B-D. The initial 
gas temperature is equal to 7000 K. 

In all models we assume that supernovae explode ran- 
domly in the disk and supply cosmic rays with energy 
equal to 10% of 10 51 erg SN kinetic energy output. We 
assume that the cosmic ray energy is injected instanta- 
neously into the ISM, as a consequence of each SN event, 
with a Gaussian radial profile (r$N = 50 pc) around the 
explosion center. Simulations of all four models exam- 
ined in this paper have been performed in a Cartesian 
domain with parameters summarized in Table [TJ 

The amplification of the regular magnetic field is iden- 
tified with the amplification of the azimuthal magnetic 
flux component averaged over all XZ slices through 
the discretized computational domain. In Figure [T] we 
present the time evolution of the azimuthal component of 
magnetic flux and the total magnetic ene rgy. Analogous 
lots for run A have been presented in lHanasz et al. I 
2004) . It is apparent in Figure Q] that the amplifica- 
tion of magnetic flux and magnetic ener gy depends on 
the re sistivity as it has been reported by lHanasz et al.~1 
( 2006). Strictly speaking the efficiency of magnetic field 
amplification grows with increasing the resistivity, within 
the assumed range of variations of the resistivity param- 
eter. Both magnetic flux and energy start to grow at 
t ~ 100 Myr and continue until a maximum (which is re- 
sistivity dependent) is reached. After reaching the max- 
imum the curves of magnetic flux and magnetic energy 
behave rather chaotically for all values of resistivity. 

In order to measure the amplification of the mean 
magnetic field we perform an averaging of B x and 
B v magnetic field co mponents as it has been done by 
lHanasz et al~l (j2004h for run A. We plot the mean mag- 
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TABLE 1 

Parameters of models examined in this paper. (*) The conversion rate, presented as 10% in[Hanasz et al. (20Q||), was in 

FACT EQUAL TO 100%, DUE TO A TRIVIAL CALCULATION MISTAKE. THEREFORE, THE OVERALL INJECTION RATE OF COSMIC RAY ENERGY WAS 
EQUIVALENT TO A REALISTIC ONE, CORRESPONDING TO SN RATE— 20 KPC~ 2 MyR _1 , WITH THE ENERGY CONVERSION FACTOR = 10%. 



Model 


A 


B | C J D 


Domain sizes [kpc] 
Resolution 

Vertical gravity at R [kpc] = 

Gas column density [cm -2 ] 

Angular velocity H [Myr -1 ] 

SN rate [kpc -2 Myr -1 ] 

Initial a = e ma g/e g as 

Diffusion coefficients Ku, K± [cm 2 s -1 ] 

Conversion rate of SN kinetic to CR energy 


0.5 x 1 x 1.2 
50 x 100 x 120 
8.5 

0.05 
2 

lO" 8 
3 x 10 27 , 3 x 10 26 
10%W 


0.5 x 1 x 4 
50 x 100 x 400 
5 

27xl0 20 
0.05 
130 
10~ 4 
3 x 10 27 , 3 x 10 26 
10% 


Resistivity r\ [cm 2 s _1 ] 


3 x 10 24 


x 10 24 | 3 x 10 24 | 30 x 10 24 




Fig. 1. — Time evolution of the magnetic flux (left plot) and magnetic energy (right plot) for the simulations B, C and D. We applied 
the normalization of both quantities to their values at t = 0. 



netic field B(z) for runs B, C and D in Figure [5] at 
different times. We find that in the growth phases of 
both magnetic flux and magnetic energy, shown in Fig- 
ure [T] (100 Myr< t <500 Myr for runs B and C, and 
100 Myr< t <800 Myr and t >1050 Myr for run D), the 
averaged magnetic field components B x (z) and B y (z) are 
apparently more regular, while at the decay phase the 
vertical fluctuations of the mean field dominate. 

A close inspection of the mean, horizontally averaged 
azimuthal magnetic field component along the simulation 
time for run B, shows that the magnitude of the mean 
field does not grow in the period Myr< t <500 Myr, 
but only the width of the central layer, of coherent (or 
uni-directional) azimuthal field, increases slightly. In the 
opposite case of run D, one can note the growth of the 
magnitude of the mean azimuthal field, till t = 800 Myr, 
as well as the growth of the width of the central layer of 
coherent azimuthal field. The radial component of the 
averaged field displays much more temporal and spatial 
fluctuations than the azimuthal component in runs B-D. 

Generally speaking fluctuations of the averaged field 
appear to be much more pronounced in runs B-D than 



in run A. This may be related to essentially more vig- 
orous dynamics, related to a significantly larger SN rate 
(at the galactocentric radius 5 kpc, see Table 1) of the 
system in simulations B-D as compared to simulation A 
corresponding to the SN rate at the galactocentric radius 
of Sun. This difference gives rise to a vertical wind, pow- 
ered by CR pressure gradient, reaching velocities of the 
order of 100 km/s at the altitude of 2 kpc. 

The oscillations present in the averaged field migrate 
out of the disk midplane at large altitudes. On the other 
hand the behavior of the averaged field is almost station- 
ary (coherent in space and time) near the disk midplane 
in the large resistivity run D. We shall associate the men- 
tioned coherence with the presence of the mean field. In 
this sense the mean field is clearly generated in runs A 
and D, while its generation in runs B and C is ques- 
tionable. The net azimuthal magnetic flux, apparent in 
the growth phases of magnetic field for run D, shown in 
Figure [TJ is a clear signature of the growth of the mean 
field. 

The set of models B-D forms an interesting sample of 
experiments to be examined from the point of view of 
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(Run B) t= 200[Myr] 

2000 I ' ' 





(Run B) t= 400[Uyr] 
2000 I ' ' 
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(Run B) t= 800 [Myr] 
20CC 1 




(Run C) t= 200[Myr] 

20001 ' ' 




(Run C) t= 400[Uyr] 
20001 ' ~" — 



(Run D) t= 600 [ Uyr ] 
2000 I r " 




(Run C) t= 600[Myr] 
20001 1 7T7 1 

<B_> ^ 




O.OD 

2<B„>, <B,> [/iG] 




(Run C) t= 800 [Myr] 
20001 1 - ~ "" ""' 



2<B,>, <B r > [ftC] 



°) t= 800[Uyr] 
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D ) t= 1000[Myr] 






2000 
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<B,> 






<B,> 




2<B,> 






2<B,> 








1000 
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.5 










-1000 










-2000 







0.00 0.26 
>, <B r > [mG] 




-0.08 0.00 0.08 

2<B„>, <B,> 1/iG] 

(Run D) t= 1100 [Myr] 
2000 I 1 - ' 




FlG. 2. — Evolution of the horizontally averaged magnetic field (solid line — azimuthal component, dotted line - radial component) for runs 
B, C and D. It is noticeable that in the growth phases of the both magnetic flux and magnetic energy, shown in Fig. ^(100 Myr< t <500 Myr 
for runs B and C, and 100 Myr< t <800 Myr and t >1050 Myr for run D), the averaged magnetic field components B x (z) and B y (z) are 
apparently more regular, while at the decay phases vertical fluctuations of the mean field dominate. 
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dynamo theories. In the forthcoming parts of the paper 
we are going to address the question to which extent the 
available dynamo theories are capable of reproducing the 
magnetic field evolution relying on dynamo coefficients. 
It is quite obvious that the truncation of the power series 
of the electromotive force, at linear or quadratic terms 
should lead to some discrepancy between magnetic fields 
measured in experiments and the ones reconstructed on 
the base of dynamo coefficients. This discrepancy can 
be directly measured and quantified in a statistic way. 
Another simple criterion for the applicability of dynamo 
theories can be proposed, however, in a more qualitative 
fashion. We shall propose a criterion, which relies on 
a simple observation wheather the growth and decay of 
reconstructed magnetic field correlates with the growth 
and decay of the original field. 

3. MEAN FIELD DYNAMO THEORY 

The evolution of magnetic fields in a dynamic system 
can be investigated in the mean field approach, where the 
separation of the total magnetic field into two parts, the 
mean large-scale and fluctuating small-scale components, 
is the fundamental assumption (|Parker I 1955; M offatt I 
119781: lRadle71ll980f lV In the mean field dynamo theory 
the main interest is to determine the influence of the 
small-scale fluctuations upon the mean magnetic field. 
Following the basic steps of the mean field theory one 
can derive the dynamo equation. 

We start from the equation of induction with the dif- 
fusion term: 

3B 

- = Vx(uxB-t|VxB). (1) 

Next, we separate the velocity and magnetic fields into 
the mean and fluctuating components: 

B=B+b, 6 = 0, (2) 

u = V + v, v=0, (3) 

where overline designates the mean field obtained from 
the averaging procedure. To obtain the mean field com- 
ponent properly, the averaging procedure should fulfill 
the required conditions called the Reynolds rules (see 
iBrandenburg fc Subramanian 1120051 e.g.): 

C7i + U» = Ui+U», U = U, Uu = Q, ZhU~2 = U±U 2 , 

(4) 

dU/dt = dU/dt, dU/dxi = dU/dxt . (5) 

After the averaging operation we get two equations de- 
scribing the evolution of the mean and the fluctuating 
parts of the magnetic field, B and b respectively: 

OB 

— = V x \vxb + VxB-rjV xS\ , (6) 

— = V x [7x6 +v xB + (v x 6 — v x b) — rjV x b] , 
at J 

(7) 

where the large-scale part of the vector product of fluc- 
tuating parts of the velocity and magnetic fields is the 
electromotive force: 

£=vxb. (8) 

The investigation of £ is the crucial point in the dy- 
namo theory. It is considered a functional of v, v and 
B, which can form nonlinear higher order moments (see 
iRadler I [1980; Brandenbur g fc Subramanian 1 120051 and 
the discussion in i j3.ll and £|3.2[) . 



3.1. Conditions for the linear approximation of the 
electromotive force 

The simplest approximation of the dynamo equations 
in order to obtain the electromotive force prescribed in 
the form of Eq. © contains the linearization of the equa- 
tions for the fluctuating quantities. The first step of ap- 
proximation is an assumption that the fluctuations are 
small, so we can neglect all nonlinear terms, which con- 
sist of fluctuating fields of orders higher than one. At 
first, this is done by neglecting the termv x b in Eq. ([7]). 
This can be true for example if R m is small, what cannot 
be fullfiled in the case of ISM. The term V x b is usually 
neglected because it describes an advection of the small 
scale magnetic field. However, this assumption cannot 
be justified in the case of the presence of a shear in the 
system. Strong shears could le ad to a new dynamo effect, 
the s hear-current effect (see iRogachevskii fc Kleeorin I 
2003]). In the case of a small R m one can neglect both 
the nonlinear terms Vx(vxb—vxb). In most astro- 
physical applications, these terms can still be neglected 
if the correlation time T corr of the turbulence is small, 
T corr v rms kf <c 1, where v rms and kf are typical veloci- 
ties and wave nu mbers associated with the random ve - 
locity field v (see IBrandenburg fc Subramanian I l2005h . 
Finally, we result in the simplified equation describing 
the evolution of the fluctuating fields: 

8b 

— = VxvxB. (9) 

dt K ' 

The absence of the diffusive term in Eqs. ©-dH]) follows 
from the assumption, that in the correlation time the 
advection strongly dominates over diffusion effects, so 
Rm 3> 0. This is valid for most of astrophysical applica- 
tions where R m ~ 10 10 — 10 20 . To calculate £ we inte- 
grate Eq. §§§ to get 6, take the cross product with v and 
average. If the mean magnetic field varies only weakly 
in space and time, £ can be represented as a function of 
B and its derivatives: 

£i = 6B - V x B = oyBj + Ajk--j (10) 

The quantities £ and B with its derivatives can be cal- 
culated directly from the simulation data. The equa- 
tion (|10p involves two tensor quantities, a = {ay} and 

$ = {Ajk}j describing the efficiencies of the amplifica- 
tion and diffusion of the mean magnetic field, respec- 
tively. Both coefficients are necessary for the mean field 
dynamo. Dynamo coeffici ents a\\ and ffiik ca n be fitted 
using statistical methods (jKowal et al. 1120051 ) . We use a 
very popular and flexible method of the multid imensional 
linear regression (see, e.g. JPress et al. 1119971 §15.4). 

3.2. Nonlinear approaches to the dynamo 
3.2.1. Blackman- Field approach 

iBlackman fc Field I (|2002f ) instead of imposing the 
form of the £, solved the dynamo equation using the 
time evolution of £ 

d t £ = d t v x b +v x dfi. (11) 

Using equations for the evolution of the small-scale 
velocity and magnetic fields (see Eqs. 4 and 5 in 
IBlackman fc Field Il2002t) and assuming the incompress- 
ibility of the flow and isotropy of the resulting velocity 
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Fig. 3. — Power spectra of the velocity (left plot) and magnetic (right plot) energy along each direction (thin solid, dotted and dashed 
lines for X-, Y- and Z-direction respectively) for run B at time 1000 Myr. The thick solid line corresponds to the Kolmogorov power law 
E ~ fc" 5 / 3 . 



and magnetic field correlations for terms linear with B, 
they obtained an equation for the time evolution of the 
electromotive force £ 

af _ 

~ f3V xB + vV 2 vxb+ Aux V 2 6 + T v 



at 



oB 



(12) 

where a = (1/3) [b • V x b - v ■ V x v) and = (1/3)1? 
are dynamo coefficients, v and A are the viscosity and 
magnetic diffusion coefficients respe ctively, and T v and 
T M are the triple correlations (see Blackm an fe Field"! 
2002, for description). 

The above and following equations are valid if the di- 
vergence of magnetic helicity flux vanishes. Neverthe- 
less, we omit the considerations of magnetic helicity in 
this paper. We intend to devote a separate paper for the 
cosmic-ray driven dynamo examined from the point of 
view of magnetic helicity conservation. 

Dynamo coefficients a and (3 do not change if we limit 
our consideration to the component of £ parallel to_B. In 
this case we obtain the following simpler expression 

d£n „B 2 ~B\7xB - 



dt 



a 



\B\ 



(3- 



(13) 



Here £ co rresponds to the macro physical dissipation 
terms fsee iBlackman fc Field 1120 02). The coefficients a 

and [3 are slightly different from the usual dynamo co- 
efficients, because they do not include any characteristic 
timescale like correlation or relaxation time, and thus 
have different units. 

In the approach described above the authors consid- 
ered all quadratic terms in the mean field in the applied 
electromotive force approximation. 

3.2.2. Kleeorin-Rogachevski approach 

The standard local thin disc approximation for dynamo 
equations for an axisymmetric magnetic field leads to the 
following equations for the mean radial field B r = R a b r 



and toroidal field B^ for the a£7-dynamo pr oblem (see 
iRuzmaikin et al. Ill988t IKleeorin et al. 1120031 e.g.) 



dbr = d(a(B)B^ 
dt 



dz 

dB^ 
dt 



d_ 

dz 



Va(B) 



db, 



d(V A (B)b r ) 



dz 



Db r 



d_ 

d?i 



d^ 

dz 



(14) 
(15) 



where D = Ru,R a = [r(dil/dr)h 2 /t7t][^Q!*/??t] is the dy- 
namo number (h is the disc thickness, a* is the maximum 
value of the hydrodynamic part of the a-effect), t)a(B) 
and t]b(B) are the nonlinear turbulent magnetic diffusion 
coefficients of poloidal and toroidal magnetic fields, the 
nonlinear function ~Va(B) = [i]a(B) — rjB(B)] d^hiB), 
and a(B) is the total nonlinear a-effect, which includes 
its kinetic and magnetic parts 



a(B) 



a 



(16) 



These quantities are determined by the corresponding 
helicities (x v and x m ) an d quenching functions (4> v and 

4>m) 

a(B)= X v MB)+X c (B)<P m (B), where B = \B\. 

(17) 

iRogachevskii fc Kleeorin I (|20000 derived the quenching 
functions for an anisotr opic turbulence, which are given 
by Eq. (10) and (11) in IKleeorin et al. I (|2000D . The ki- 
netic helicity is given by x v — — ( r /3)(f ■ (V xv)). The 
magnetic part of the a-effect is introduced by the func- 
tion \ C (B ) determined by the e volutionary equation (see 
Eq. 12 in lKleeorin et al. Il2000f ). 

Eq. (fT7l) with helicitie s and quenching functions de- 
rived by IKleeorin et al~l (|2002f ) contain the main non- 
linearities of a-effect. They include the conventional 
and algebraic quenching and dynamic nonlinearity de- 
sc ribed by the evolution of the magnetic helicity (Eq. 12 
in lKleeorin et al. Il2002l) . 
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Fig. 4. — Power spectra of the velocity (left plot) and magnetic (right plot) energy along each direction (thin solid, dotted and dashed 
lines for X-, Y- and Z-direction respectively) for run D at time 1000 Myr. The thick solid line corresponds to the Kolmogorov power law 
E ~ fc- 5 / 3 . 
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Fig. 6. — Time variances of the mean values of ratios V X (p X 
6)/V X (u xB) (solid line), V X (v x6)/V X (v xB) (dashed line) and 
V X (V"xf»)/V X (p xB) (dash-dot line) for experiments A (upper 
plot) and B (lower plot). A condition for a linear approximation 
of £ requires these ratios to be much smaller than 1, however this 
condition is not fulfilled for the presented experiments. 

4. RESULTS 

4.1. Tests of the fulfillment of the basic assumptions of 
dynamo theory 

4.1.1. Averaging method and verification of Reynolds rules 

The choice of method for the field decomposition into 
its large and small-scale parts is a very important step 
in the practical application of the mean field dynamo 
theory. The proper averaging procedure should fulfill 
the Reynolds rules described by Eqs. (JI|)-([5|). Otherwise, 
the assumptions under which we separate the induction 
equation into two equations describing the evolution of 
the mean and fluctuating parts of-B could be violated. 

Due to the different boundary conditions imposed at 
different boundaries of our simulation domain (periodic 
at the azimuthal boundaries, periodic with shear at the 
radial boundaries and open at the vertical boundaries), 
we use the simplest averaging method: the averaging over 
the horizontal planes, as it has been done by Hanasz et 
al. 2004). The other more advanced methods, e.g. Gaus- 
sian smoothing, would require dividing the procedure of 
averaging into two steps: the backword shearing of the 
domain to get back a full periodicity at the azimuthal 
boundaries and the actual averaging. In the case of av- 
eraging over the horizontal planes, the former step can 
be omitted. 

4.1.2. Verification of scale separation and anisotropy in 
turbulent spectra 

One of the basic assu mptions of the l inear and non- 
linear dynam o theories (|Parker I 119551 : iMoffatt I 119781 : 
iRadler |[l980l . see also is the scale separation, which 
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Fig. 7. — Time evolution of the kinetic turbulent (dashed line) 
and the small-scale magnetic (solid line) energies for three models: 

A, B, and C (top, middle and bottom plots respectively). 

means that all quantities like the velocity and mag- 
netic fields have two characteristic time and length 
scales: the large and the small ones. Such scales 
should be apparent as separate peaks in the velocity and 
magnetic spectra, provided that the computational do- 
main is large enough, showing the characteristic scale 
lengths of waves in a turbulent reg ion (see Fig. 8.1 in 
iBrandenburg fc Subramanian 1120050 . 

In Figured] we show power spectra of kinetic and mag- 
netic energy fluctuations taken at t =1000 Myr for model 

B. We show only one time snapshot because in the case of 
vanishing explicit resistivity (n = 0) the power spectra of 
magnetic energy are very similar all over the simulation 
time. The kinetic energy spectra are obtained on the ba- 
sis of density and velocity component distributions after 
subtracting the large scale sheared azimuthal fluid flow. 
In order to compare the spectra in our models to the 
Kolmogorov spectrum we present also a line represent- 
ing the k -5 / 3 slope in the plots. The presented spectra 
of kinetic energy exhibit in the large wavenumber limit a 
similarity to the Kolmogorov k -5 / 3 spectrum. The mag- 
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netic spectrum appears to be much flatter in the same 
wavenumber range. It is apparent that both the kinetic 
and magnetic spectra are monotonic, thus we can state 
that there is no separation in the turbulent scales of ve- 
locity and magnetic field components. 

The fact that the spectra show the dominant power 
on the box scale has plausibility for the following rea- 
sons: (i) The models rely on the buoyancy of cosmic 
rays and magnetic fields. Cosmic rays, even if injected 
locally, diffuse effectively along the mean magnetic field 
to fill the box scales quickly, (ii) Magnetic loops recon- 
nect due to an explicit or numerical resistivity. As it has 
been originally suggested by Parker (1992), the discussed 
cosmic-ray driven dynamo relies on the inverse turbulent 
cascade. The box sizes, signifying the largest available 
spatial scales, are limited by the simulation times, which 
are equal to several months of CPU time in single pro- 
cessor runs for all presented simulations. 

In Figure 0] we present analogous spectra for model D 
with large resistivity rj =10 pc 2 /Myr, at t = 1000 Myr. 
We can see that the experiment with larger resistivity 
gives similar kinetic energy spectrum and a steeper mag- 
netic spectrum, which are both close to k -5 / 3 spectrum, 
than model B. This means that the magnetic energy re- 
sides more at larger scales in models with larger resistiv- 
ity, according to expectations. 

The differences in spectra of magnetic field and veloc- 
ity components can be discussed in terms of correlation 
length of these components along x, y and z axes. The 
temporal evolution of correlation length of all the men- 
tioned components in each direction, for model B, is dis- 
played in Figure [5J It is apparent in Figure [5] that the 
correlation lengths of v x and v y are in all three directions 
close to 100 pc, while the correlation length of v z is about 
70 pc in x and y direction and oscillates around 350 pc 
in z direction. The large correlation length of v z along 
z direction is a manifestation of a vertical wind powered 
by the vertical cosmic ray pressure gradient. In the case 
of magnetic field the correlation length of B x in all di- 
rections is comparable to 50 pc, the correlation length of 
By varies around 40 pc in the x direction, around 70 pc 
in the z direction, and varies from 100 pc to more than 
400 pc, and again down to 100 pc during the system 
evolution. Similarly to v z the correlation length of B z 
along the z direction is much larger than the correlation 
length of B z in the x and y directions. This is again a 
signature of vertical, nonuniform wind, which stretches 
the magnetic field in the vertical direction. 

4.1.3. Tests of validity of the linear approximation 

For the next test we check the relative magnitudes of 
all terms in the £ calculated from our experiments A and 
B (see £J3J) - The ratios of electromotive force terms cal- 
culated according to Eq. are presented in Figure [6] 
for models A (top) and B (bottom). As it was shown 
in £13.11 the linear dynamo theory neglects all quadratic 
terms of the mean field in the EMF_ equation so they 
remain only with the term V x (v x B) (sec Eq. |9]). To 
check how large are the three other terms (V x (» x 6), 
V x (vxb, V x (Vxb)) in comparison with this term we 
divide them by the first one. It is apparent that after 
10 Myr for the model A the highest ratio is obtained for 
the term V X (v X b) , which grows even ten times higher 
than the term V x (»xB). The second term V x (Vxb) 
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Fig. 8. — Comparison of the time evolution of the component 
of £ parallel to the local mean magnetic field reconstructed by the 
Blackman approach for three values of f (0.0 — dotted line, 0.01 — 
dashed line, 0.05 - dash-dotted line) to £m taken directly from the 
experiments (solid line) for model A for Z=400pc. 

grows few times higher than the first term. The third 
term V x (v x b) shows ratios changing from 0.1 to a few 
times larger. The evolution of term ratios in model B 
(Fig. [6j bottom plot) is very similar to the evolution of 
two other models C and D, so we present only this case 
in the Figure. In Figure [6] (bottom) we can see that 
in the beginning all terms studied are close to the term 

V x (v xB). Later on, the term V x (vxb) has similar 
value to V x (v xB) . The ratio of the terms V x (u xft) and 

V x (Vxb) changes their magnitude from about 0.8 in the 
beginning to above 10 in the end of the evolution. Such 
high values of the normally neglected EMF terms mean 
that in our numerical experiments of the Parker insta- 
bility, driven by the cosmic rays and shear, the resulting 
electromotive force is not in a linear regime. For these 
reasons one should consider inclusion of all the terms 
in the calcul ations of the dynamo coef ficients, as it was 
proposed by iBlackman fc Field I (|2002f ) in their calcula- 
tions of the quenched evolution of EMF. Although as it 
will appear later, the quality of the approximation of the 
mean electromotive force will not become satisfactory. 

The next issue is related to the assumption that the ki- 
netic turbulent energy is much la rger than the energy o f 
the small-scale magnetic field (see lKleeorin et al. 1120031) . 
Figure [7] presents that in model A, since 1000 Myr, both 
energies are comparable. In the case of experiments B 
and C the situation is different, showing that the mag- 
netic turbulent energy is smaller than the kinetic one 
during the whole period of calculations. For this rea- 
so n in calculations of the dynamo coefficients according 
to iKleeorin et al. I (|2003h we apply their approximations 
for models B-D. 

4.2. The fitting of the different dynamo coefficients to 
the calculated EMF 

4.2.1. Implementation of the Blackman- Field approach 

The next goal is to fit the dynamo coeffic i ents c alcu- 
lated with the help of the IBlackman fc Field I (|2002f ) pre- 
scription to the electromotive force r esulting from our nu - 
merical models A, B and C (see also lKowal et al. 1 12005') . 
The authors (BF02) applied the nonlinear dynamo the- 
ory including all quadratic terms in the mean field in 
the EMF approximation (see i )3.2.1|) . We shall check if 
their approximation gives similar evolution of the EMF 
to the original one. We apply Eq. (fT"3|) to obtain the time 
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Fig. 9. — Comparison of the time evolution of the component 
of S parallel to the local mean magnetic field reconstructed by the 
Blackman approach for three values of £ (0.0 - dotted line, 0.01 - 
dashed line, 0.05 — dash-dotted line) to £j| taken directly from the 
experiments (solid line) for models B and C (upper and lower plot, 
respectively) for Z=400pc. 

derivative of the component of the EMF parallel to the 
local magnetic field, the values of coefficients a and (3 
and the remaining terms in that equation. Figs. [5] and 
[5] present the time evolutions of the original and the re- 
constructed EMF for models A, B and C in the chosen 
time periods. 

The large discrepancies of the original and recon- 
structed EMF signify a weak application of the Black- 
man model to the results obtained in the simulations. 
We shall further discuss in H4.3.1I the results of recon- 
struction of mean magnetic field components, based on 
Blackman k Field approach. 

4.2.2. Implementation of the Kleeorin-Rogachevski approach 

The second analysis concerns the pre scription of the 
dynam o coefficients given in the paper iKleeorin et al. I 
(|2003f ) (see M3.2.2|) . The authors apply in the electro- 
motive force prescriptions only the term v xb and ne- 
glect the terms V x b, v x B, and v x b, but they 
take into account the anisotropy of turbulence (s ee also 
iRogachevskii k Kleeorin~ll200lD . IKleeorin et "all (|2003f ) 
applied their quenched form of the dynamo coefficients to 
the large-scale dynamo simulations. In £|3.2.2l we present 
the set of equations from their paper, which we include 
in our calculations. Again we use our model B in order 
to calculate both magnetic and kinetic a, as well as the 
different forms of the diffusion coefficient. We add all re- 
sulting terms to g et the componen t s of th e electromotive 
force according to IKleeorin et al~l (|2003l ). 

The time evolution of the actual (the solid line) and 
reconstructed (the dashed line) electromotive force com- 
ponents averaged over the planes are presented in Fig- 
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Fig. 10. — Comparison of the time evolution of X-, Y-, and 
Z-component of £ (top, middle and bottom plots respectively) re- 
constructed by the Kleeorin-Rogachevski approach (dashed lines) 
to those taken directly from the experiments (solid lines) for ex- 
ample model B for Z=400pc. 

ure [10] One can see that differences between the ac- 
tual and the reconstructed electromotive forces are re- 
markable, however, their time averaged values seem to 
match the reconstructed electromotive forces, at least in 
the y direction. The discussion of reconstructed mean 
magnetic field basing on the Kleeorin k Rogachevski ap- 
proach will be presented in £14.3.21 

4.3. Integration of the mean magnetic field from EMF 

In the previous section we show that, the Blackman- 
Field approach does not reconstruct the electromotive 
force properly, while the Kleeorin-Rogachevski approach 
provides in some cases electromotive forces that follow 
approximately time-averaged electromotive forces taken 
from experiments. Nevertheless, these approaches per- 
mit for obtaining some useful informations about the 
electromotive force and its properties. For instance, in- 
tegrating Eq. ^ we can obtain a qualitative information 
about the amplification of the mean magnetic field. We 
also can analyze, how much of the mean magnetic field 
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can be restored by the integration of the electromotive 
forces reconstructed by each approximation. 

We perform the time integration of Eq. © substituting 
the term v x b by either original or reconstructed electro- 
motive force £ : 

an 

— = V x [£ + VxB-tjV xB] . (18) 

In addition, we could neglect the second term incorpo- 
rating the large-scale fields V and B, because in general, 
these field should be weakly varying in space if the pro- 
cedure of separation undergo the Reynolds rules. How- 
ever, in our analysis, we incorporated the plane averaging 
procedure, which could give substantially large vertical 
derivatives, thus in the above equation we take the large- 
scale component of velocity V directly from models. In 
this way, the only unknown variable is the large-scale 
component of magnetic field B. If assumptions of the 
mean field theory are fullfiled, the Eq. (fTS"]) should recon- 
struct the mean magnetic field correctly. The numerical 
integration is performed with t he use of the fourt h order 
Runge-Kutta method (see e.g. iPress et al. Ill997t ). 

4.3.1. Reconstruction from the Blackman-Field approach 

In the Blackman-Field approach the information about 
the electromotive force is incomplete, because this ap- 
proach reconstructs only the component of the £ parallel 
to the local mean magnetic field B. The perpendicular 
component of £ is lost and cannot be recovered in this 
approach. Moreover, the £\\ is now a scalar field. Due 
to this limitation we need to add an additional step in 
which all three X, Y and Z-components of £\\ must be 
reconstructed. The simplest way to perform a projection 
of the £\\ on the direction of B taken directly from the 
simulations. We should note, that by taking only £u into 
integration, even that calculated directly from the sim- 
ulations, does not guarantee a perfect reconstruction of 
theB. 

The resulting mean magnetic fields integrated from 
EMFs, which were reconstructed by the Blackman-Field 
approach are shown in Figure 1111 This figure presents 
the integrated mean magnetic field for all four models: 
A, B, C, and D (from top to bottom) for three different 
coordinates Z=0, 150 and 400pc (left, middle, and right 
columns, respectively). We see that all integrations give 
different evolutions, which are not always in the agree- 
ment with the original B. For instance, in Figure [TT] we 
see, that the integration of the EMF taken directly from 
models A-D delivers a good agreement with the original 
B only during the short initial period of the evolution. 
Later on, the reconstruction of £u from the Blackman- 
Field approach provides significantly larger growth of the 
mean fields, especially for z—0 and 150pc. When the al- 
titude is higher, the dashed lines corresponding to the in- 
tegration of the Blackman-Field's £\\ show much larger 
discrepancies than the solid lines for integration of the 
original £u , but they still unveil the growth of the mean 
magnetic field B. In models C and D (two lower rows in 
Figure ITT]) . the better compatibility with the original B 
comes from the integration of £ u taken from the simula- 
tion. In these two models, the Blackman-Field approach 
simply overamplifies the authentic evolution of the mean 
magnetic field, producing up to two orders of magnitude 



larger B. For higher altitudes (plots on the right col- 
umn), the evolutions are again comparable, although we 
still see a clear tendency to overamplify the magnetic 
field by the Blackman-Field approach. The model B is 
exceptional, because the amplification of B takes place 
only in the initial stage of the evolution. Later on, the 
strength of the mean magnetic field declines. Despite 
that fact, the integration of the actual £ gives almost 
constant mean magnetic field at the midplane, while only 
a small growth of£? at the initial stage at higher heights. 
From time T=400Myr we do not see any essential ampli- 
fication of-B. The Blackman-Field approach gives again 
the amplification of magnetic field at the initial stage. 
Later on, it varies or even declines, what is well seen at 
higher altitudes. 

The above analysis draws a conclusion that the 
Blackman-Field approach reconstructs the £ || preserving 
its property of amplification of the magnetic field. How- 
ever, this approach often overamplify the mean magnetic 
field, even in the situations when the strength of the ac- 
tual magnetic field decreases. Moreover, we see another 
conclusion related maybe not directly to the Blackman- 
Field approach, but rather to the properties of the elec- 
tromotive force. Model A shows that the parallel com- 
ponent of the electromotive force contributes to the am- 
plification of_B later than the actual amplification is ob- 
served. This is justified by the fact, that the electromo- 
tive force reconstructed from its component parallel to 
B is incomplete. This could indicate, that the amplifica- 
tion of B during the initial stage comes mainly from the 
remaining perpendicular component of £ . 

4.3.2. Reconstruction from the Kleeorin-Rogachevski 
approach 

The mean magnetic fields integrated from the electro- 
motive force reconstructed by the Kleeorin-Rogachevski 
approach are shown in Figure [12] for the same models and 
altitudes as previously. Here, the result of integration is 
more consistent with the actual evolution of-B, because 
we integrate full £ , and not only its component parallel to 
B. In the top row of Figure [T2l we see the mean magnetic 
fields integrated for model A, which we also take into 
account in order to be consistent with the BF02 presen- 
tation. The Kleeorin-Rogachevski approach reproduces 
the mean magnetic field relatively well at z=0pc, how- 
ever, only for the initial several hundred Myrs of evolu- 
tion. Later on, the discrepancies are larger. The growth 
rate of B resulting from this approach starts to decrease, 
while the growth rate of the actual field is relatively con- 
stant. At the altitudes 150 and 400pc the growth ofB is 
very consistent with that integrated from £ up to time 
about 1500pc. Later, the approach gives again weaker 
magnetic field, by about one orders of magnitude. 

In model B (second row in Fig. [T2"|) , the mean mag- 
netic field integrated from the original £ and the mean 
magnetic field taken directly from simulations are more 
or less in agreement in the whole time period, while the 
reconstructed curves show much larger discrepancies, es- 
pecially at z=0pc. This model is less interesting, because 
we do not observe an amplification of B here. In model 
C the reconstruction of the mean magnetic field given 
by the original £ is really very good, while the curves 
reconstructed according Kleeorin-Rogachevski approach 
are good only for the height z=0pc. At two higher z (150 
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Fig. 11. — Time evolution of the absolute values of the mean magnetic field integrated from original (solid lines) and reconstructed 
(dashed lines) £ for three vertical positions Z=0, 150, 400pc (left, middle and right columns, respectively) for models A, B, C and D (from 
top to bottom row). The reconstructed £ was calculated using the Blackmann-Field approach. We also show the original mean magnetic 
field taken directly from simulations (dotted lines). 



and 400pc) show that this approach gives much lower 
values of the mean magnetic field. This approach fails 
completely also in model D, where the large resistivity re- 
sponsible for the strong diffusion leads to a fast growth of 
B. Here, the Kleeorin-Rogachevski reconstruction gives 
only marginal amplification of the mean magnetic field. 

5. DISCUSSION 

In the present project the model of the CR-driven 
galactic dynamo is analyzed in order to determine its lin- 
ear or non-linear character. We analyze the conditions 
for validity of linear approximation of the dynamo equa- 
tion in our model. First, we check the main assumptions 
of the dynamo theory: the Reynolds rules, the separation 
of the large and small scales - both magnetic and kinetic. 
Then we compare the ratios of different terms appearing 
in the equation for time evolution of the fluctuating part 
of the magnetic field. We applied averaging over the hor- 
izontal planes. Further investigations revealed that the 
condition of scale separation is violated in our numer- 
ical model. In this sense our models do not fulfill the 
assumptions of dynamo theory (see As far as the 



scale separation is concerned, our magnetic and kinetic 
spectra computed for all chosen numerical models (A, 
B, C and D) are flat and are characterized by the same 
slope at all scales. This means that in our experiment no 
scale separation occurs. The next essential problem con- 
cerns the comparison of the different EMF terms. The 
plots shown in ij4.ll indicate the similarity among all ex- 
periments. The terms, which are normally neglected in 
the linear dynamo theory should be taken into account 
(see BF02 for comparison). For this reason the discrep- 
ancy between actual EMF and the reconstructed EMF is 
present in majority of experiments. 

We have check subsequently, how the dynamo non- 
line ar approximations pre sented in the literature (BF02 
and lKleeorin et aT~ll2003f ) fit the electromotive force ob- 
tained from our models (A, B, C and D). We compute 
BF02 electromotive force and present that it does not 
follow the time-evolution of the EMF obtained directly 
from our models. In particular, their values diverge from 
the original ones, and ensure too fast growth of energy 
of the reconstructed mean magnetic field. The discrep- 
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Fig. 12. — Time evolution of the absolute values of the mean magnetic field integrated from original (solid lines) and reconstructed 
(dashed lines) £ for three vertical positions Z=0, 150, 400pc (left, middle and right columns, respectively) for models A, B, C and D 
(from top to bottom row). The reconstructed £ was calculated using the Kleeorin-Rogachevski approach. We also show the original mean 
magnetic field taken directly from simulations (dotted lines). 



ancies may result from the truncation of the expansion 
series and a lack of the scale separation, the anisotropy 
of the turbulence (see §4.ip . or the compressibility of the 

S as - . . 

The approxima tions of IRogachevskii fc Kleeorin I 

(120001 120011 120031) a re used by many author s (e.g . 
Kleeorin et al. I 120031: iBrandenburg fc Sandinl 120041 
Brandenburg et al. T I2005D . That is why we also check 
if their anisotropic and quenched dynamo co efficients 
fit ou r modeled EMF. The authors applied iRadler I 
(1980) prescriptions, but neglected three of four terms 
of the EMF. It is easy to notice (see §4.2[) that the 
electromotive force obtained from our numerical model 
of the cosmic ray driven dynamo is different from the 
Kleeorin one. The obtained growth of the energy of the 
mean magnetic field reconstructed from the Kleeorin- 
Rogachevski approach is satisfactory only in models 
A and B and is too slow in other models with larger 
resistivity in comparison with numerical experiments. 

As we mentionned already, we did not discuss the ef- 
fects of magnetic hclicity transport and its relation to the 



alpha effects, although the terms related to the magnetic 
helicity current have been included in both Blackman- 
Field and in Kleeorin-Rogachevskii approaches through 
the magnetic part of the a-tensor. Considerations of 
magnetic helicity conservation is planned as an extension 
of the present work. The role of magnetic helicity cur- 
rent in o ur dynamo model can b e anticipated from the 
paper by IShukurov et al~~l (|2006f) . where helicity losses 
through galactic fountain flow are examined in context 
of catastrophic quenching of galactic dynamos. These 
authors demonstrate that a vertical galactic wind can 
advect magnetic helicity and make it possible for the 
large-scale magnetic field to reach the strength of 10 % 
of the equipartition value. 

The simmilarity of our model and the galactic foun- 
tain model is quite obvious. In both models convective 
motions are present, although their origin is different. In 
our case cosmic ray buoyancy, and in fountain model the 
buoyancy of hot gas drive convective motions reaching 
vertical speeds of 200 km/s. We point out to the appar- 
ent advection of magnetic fluctuations at large altitudes, 
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which indicate the possibility of carrying the magnetic 
helicity current. Furthermore, the qual itative behaviour 
of our models is simmilar to those by IShukurov et al.1 
(2006). The characteristic growth and decay phases 
are apparent in our dynamo models with low resistivity 
and in galactic fountain models with low vertical winds. 
The saturation of magnetic field growth at about 10% 
of equipartition value is also a feature of both types of 
models for larger resistivity values and moderate verical 
winds respectively. 

6. CONCLUSIONS 

Our paper presents the analysis of the nonlinear elec- 
tromotive forces in the model of the cosmic-ray driven 
dynamo in the galactic disk. Our results may be sum- 
marized by the following conclusions: 

1. Neither the velocity nor the magnetic field scale 
separation occurs in our model. 

2. The electromotive forces in the cosmic-ray driven 
dynamo model are nonlinear, but none of the two 



examined nonlinear approaches is capable of repro- 
ducing electromotive forces in the numerical exper- 
iments correctly. 

3. Various nonlinear prescriptions of the dynamo coef- 
ficients have been proposed by other outhors, how- 
ever, they are not capable of reconstructing the 
electromotive force resulting from experiments of 
cosmic-ray driven dynamo. Moreover the recon- 
structions of the magnetic field produce too fast or 
too slow growth of the magnetic energy in compar- 
ison with the results of cosmic-ray driven dynamo 
numerical experiments. 

Extension of the present work, including considerations 
of magnetic helicity conservation will be presented in the 
forthcoming paper. 

This work was partly supported by the Polish Commit- 
tee for Scientific Research from the grants 1 P03D 004 
26 and 1 P03D 002 28. 
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